Characterization of turbulence in inhomogeneous anisotropic flows 
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We introduce a global quantity 8 that characterizes turbulent fluctuations in inhomogeneous 
anisotropic flows. This time-dependent quantity is based on spatial averages of global velocity 
fields rather than classical temporal averages of local velocities. S(t) provides a useful quantitative 
characterization of any turbulent flow through generally only two parameters, its time average 8 
and its variance 82. Properties of 8 and 82 are experimentally studied in the typical case of the von 
Karman flow and used to characterize the scale by scale energy budget as a function of the forcing 
mode as well as the transition between two flow topologies. 

PACS numbers: 47.27.-i Turbulent flows, 47.27.Ak Turbulent flows-Fundamentals, 47.27.Jv High- Reynolds- 
number turbulence, 47.20.Ky Nonlinearity, bifurcation, and symmetry breaking 
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I. INTRODUCTION 

The classical theory of homogeneous isotropic turbu- 
lence applies to flows with zero velocity average so that all 
the turbulence properties are contained in the temporal 
velocity fluctuations [l[. Real turbulent flows, with 11011 
trivial boundary conditions and external forcing, often 
possess a non-zero mean flow that may or may not be 
stationary [26]. Quantifying the influence of turbulent 
fluctuations on the flow dynamics has been the corner 
stone of turbulence theory. In early experimental works, 
velocity measurements were performed through hot wire 
probes providing temporal variations of the velocity at 
fixed points. A natural parameter quantifying the inten- 
sity of turbulence has therefore been introduced as: 



<v 2 -v 



V v 

where V is the velocity and X refers to the time average 
of variable X. While such a parameter clearly character- 
izes fluctuations in a homogeneous flow, one may ques- 
tion its relevance in more general anisotropic inhomoge- 
neous flows where turbulence intensity depends on the 
measurement point (cf. Figs. []Ja) and (d)). For exam- 
ple, i diverges around stagnation points or shears so that 
no robust global quantity can be built by integrating or 
averaging i over the whole flow. This question is receiving 
an increasing interest since the advent of Particle Image 
Velocimetry (PIV) measurements which provide instan- 
taneous snapshots of the velocity field. An illustration of 
the importance of that issue may be given in the typical 
inhomogeneous anisotropic case of the von Karman flow. 
This flow, generated in between two counter-rotating co- 
axial impellers, has received a lot of interest @, H, 0, IH, @| 
as a simple way to obtain experimentally a very large 
Reynolds number flow in a compact design (Re ~ 10 6 



in a table top water apparatus). In the equatorial shear 
layer of such a flow, fluctuations are large and exhibit 
similar local properties as in large Reynolds number ex- 
perimental facilities devoted to homogeneous turbulence 
[U 0, H, @] ■ Away from the shear layer, one observes a 
decrease of the turbulence intensity [l(J, [TTJ] as seen in 
Figs. Ola) and (d). Overall, the flow is strongly turbu- 
lent, so that the instantaneous velocity fields, measured 
by means of a PIV system, strongly differ in a non-trivial 
manner from their time average (cf. Figs. |Tfb) and (c)). 

First, we introduce in Sec. |TT]a global time-dependent 
quantity 5(t) that allows to characterize globally turbu- 
lent fluctuations in any inhomogeneous anisotropic flow. 
In Sec. Mil we present our experimental setup and show 
how to compute practically S(t) from Stereoscopic Par- 
ticule Image Velocimetry (SPIV) measurements. In Sec. 
IIV1 we first show that 6(t) provides a useful quantitative 
characterization of turbulent flows through generally only 
two parameters, its time average 5 and its variance 62, 
that quantify respectively the level of fluctuations com- 
pared to the mean flow and their ability to disturb the 
mean flow. These two quantities are introduced as a 
generalization of the classical local turbulence intensity 
i. Finally, we discuss applications of these new global 
measures of turbulence intensity: we show how these pa- 
rameters may be used to study the influence of the forcing 
conditions on the flow. 



II. GENERALIZATION OF THE TURBULENCE 
INTENSITY 



We consider a general anisotropic inhomo- 
geneous turbulent flow, with velocity field 
V(x,t) = [u(x,t),v(x,t),w(x,t)]. In order to smooth out 
the possible spatial inhomogeneities of the turbulence 
intensity, we define a scalar quantity through : 
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FIG. 1: Maps of turbulence intensity i, (a) and (d), corresponding time averaged, (b) and (e), and typical instantaneous, (c) 
and (f), velocity fields for two von Karman flows ((a-c) and (d-f)) at the same Reynolds number but with slightly different 
geometry. Both flows correspond to the same experimental setup — TM60(+) — , except for a narrow annulus inserted in the 
median plane of the flow in (d-f). In (a) and (d), for better visibility, the gray scale has been saturated at i = 6 since values 
up to i ~ I0 4 can be reached in the median plane due to weak mean velocities and large fluctuations. Additionally, vectors of 
the instantaneous velocities have been scaled by half with respect to the mean field vectors. 

turbulent flow, 8 = i 2 + 1, so that 5 is a global general- 
ization of the local turbulence intensity; second, for real 
flows such as the von Karman flow, 8 contains an addi- 
tional information regarding how far the instantaneous 
flow is from the mean flow. Indeed, 

(5-l)(V 2 ) = {( V - V) 2 ) 

where the sum runs over the volume V of the flow. The 
quantity under the overline is nothing but the square of 
the mean distance (using norm 2) between the instan- 
taneous flow and the mean flow in the functional space. 
Therefore, the quantity 5 — 1 measures how far, on av- 
erage, the instantaneous velocity field is from its time 
average. If <5 is close to one, one therefore expects the 
instantaneous flow to strongly resemble the mean flow. 
If 8 is much greater than 1, the instantaneous field will 
be more remote from the mean flow. This interpretation 
of 8 will be used in appendix lAl to draw a rough study of 
the convergence of the von Karman flow toward its time 
average. 



where (X) and X refer respectively to spatial and time 
average of X, and V(x, t) 2 = u 2 + v 2 + w 2 is the local 
kinetic energy density at time t. The quantity 8 can also 
be viewed as the ratio of the total kinetic energy of the 
instantaneous flow to the total kinetic energy of the mean 
flow: 

W E (V) W 

Note that this scalar parameter is time dependent and 
generally widely fluctuates in time (cf. Fig. [2]). We then 
define two time independent parameters, 8 and 82, that 
are respectively the time average and the variance of 8(t) 



(V 2 (t))/(V ), 



S(t) 2 -8(t) . 

We show in Sec. IHII that these two quantities fully char- 
acterize 8(t) provided that the considered time series of 
sampled fields are uncorrelatcd. Two interesting physical 
interpretations of 8 can be drawn: first, in a homogeneous 
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FIG. 2: Example of the temporal evolution of S(t) at two 
time scales: a full record of about 20 minutes on the left, and 
a focus on a 40 seconds sample on the right. A schematic 
illustration of the meaning of S and 82 is also provided over 
the longest time record. 



III. EXPERIMENTAL SETUP AND DATA 
PROCESSING 

A. The von Karman flow 

1. Experimental setup 

In order to illustrate and apply these concepts, we have 
worked with a specific axisymmctric turbulent flow: the 
von Karman flow generated by two counter-rotating im- 
pellers in a cylindrical vessel. The cylinder radius and 
height are respectively R = 100 mm and 240 mm. We 
have used two sets of impellers named TM60 and TM73 
[I?) . These two models are flat disks of respective diam- 
eter 185 and 150 mm fitted with radial blades of height 
20 mm and respective curvature radius 50 and 92.5 mm. 
The inner faces of the discs are H = 180 mm apart. 
Impellers are driven by two independent motors rotating 
up to typically 10 Hz. More details about the experimen- 
tal setup can be found in [l2j . The motor frequencies can 
be either set equal to get exact counter-rotating regime, 
or set to different values /1 7^ /2- We define two forcing 
conditions associated with the concave (resp. convex) 
face of the blades going forward, denoted in the sequel by 
sense ( — ) (resp. (+) ). We also work with two different 
vessel geometries, allowing the optional insertion of an 
annulus — thickness 5 mm, inner diameter 170 mm — in 
the equatorial plane. Both the forcing condition and the 
annulus insertion strongly influence the level of fluctua- 
tion in the flow, thereby allowing to test experimentally 
the sensitivity and the relevance of 8(t), 5 and 62- The 
working fluid is water. 



2. Control parameters 

From the two motor frequencies, /1 and /2, we define 
two control parameters: a Reynolds number, Re, and a 
rotation number, 8. For the experiments described in 
this paper, the Reynolds number, 

Re = v{h+ f 2 )R 2 v-\ 



where v is the fluid kinematic viscosity, ranges from 
1.25 x 10 5 to 5 x 10 5 so that the turbulence can be con- 
sidered fully developed. The rotation number, 



measures the relative influence of global rotation over 
a typical turbulent shear frequency. Indeed, the exact 
counter-rotating regime corresponds to 8 = and, for 
a non-zero rotation number, our experimental system is 
similar, within lateral boundaries, to an exact counter- 
rotating experiment at frequency / = (/1 + /2V2, with 
an overall global rotation at frequency — /2)/2 fa. ITol] . 
In our experiments, we vary 8 from —1 to +1, exploring 
a regime of relatively weak rotation to shear ratio. 



3. Mean flow topology 

In the exact counter-rotating regime, i.e. for 8 = 0, the 
standard mean flow is divided into two toric recirculation 
cells separated by an azimuthal shear layer (cf. Fig. QJb) 
and (e)). As 8 is driven away from zero, a change of 
topology occurs at a critical value 8 C : the mean flow bi- 
furcates from the two counter-rotating recirculation cells 
to a single cell Q. 8 C depends on the forcing and the 
geometry. For example, in the configuration with TM73 
impellers, rotation sense (+) and the annulus, we mea- 
sure 8 C = 0.17 ± 0.01 through torque measurements. In 
the case TM60(— ), this turbulent bifurcation becomes 
highly singular and gives rise to multistability between 
the two turbulent flow symmetries, 0(2)/SO(2) 



B. Measurements and data processing 

Measurements are done with a Stereoscopic Particle 
Image Velocimetry (SPIV) system. The SPIV data pro- 
vide the radial, axial and azimuthal velocity components 
on a 95x66 points grid covering a whole meridian plane 
of the flow through a time series of about 5000 regularly 
sampled values. The sampling frequency is set between 1 
and 4 Hz, corresponding to one sample record every 1 to 
10 impeller rotations. The total acquisition time is about 
ten minutes, i.e. one order of magnitude longer than the 
characteristic time of the slowest patterns of the turbu- 
lent flow. Fast scales are statistically sampled. 
The velocity fields are non-dimensionalized using a typi- 
cal velocity Vo = 2-7ri?(/i + /2V2 based on the radius of 
the cylinder and the rotation frequencies of the impellers. 
The resulting velocity fields are windowed so as to fit to 
the boundaries of the flow and remove spurious velocities 
measured in the impellers and the boundaries. The re- 
sulting fields consist of 58x58 points velocity maps. Two 
types of filtering arc further applied to clean the data: 
first, a global filter to get rid of all velocities larger than 
3 x Vb; then, a local filter (based on velocities of nearest 
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neighbours) to remove isolated spurious vectors. Typi- 
cally, 1 % of the data are changed by this processing. 



C. Data analysis 

We use two different methods to compute S(t) from our 
PIV measurements. In the direct method, we compute 
5 by spatially averaging the kinetic energy density of in- 
stantaneous and time averaged flows. Since we measure 
the full velocity field in a single meridian plane only, we 
compute 3D spatial averages of any quantity X assuming 
the statistical axisymmetry of the von Karman flow such 
as: 



(X(r,z)) = 



1 



R 2 H 



R pH/2 

r\dr I dz X(r, z). 



-R 



-H/2 



Since velocity fields are discrete in space and time, spa- 
tial integration is done with a classical numerical trape- 
zoid summation method whereas time integration is per- 
formed through simple summation. The parameter S can 
also be obtained by a spectral method. For this, we com- 
pute the spatial Power Spectral Density (PSD) of the 
discrete velocity fields as: 

E(k r ,k z ,t) = WW*, 

where W — W(k r , k z ,t) is the two dimensional spatial 
Fourier transform of W(r, z, t) = y/\r\ V(r, z,t) and W* 
its complex conjugate. Then, we compute 6 using Parsc- 
val's identity: 



dk r dk z WW* 



and get: 



5(t) 



J J dk r dk z E(k,t) 



J J dk r dk z E (k) 
where Eq is the PSD of the time averaged flow 
E (k r ,k z ,t) =WW*. 

IV. RESULTS AND APPLICATIONS 



(2) 



and geometry conditions studied in the present paper 
(cf. Fig. EDJb)). A striking feature is that global rotation 
does not change the statistical properties of the flow since 
the PDFs are Gaussian at any 9. Consequently, we can 
describe and fully characterize the total energy temporal 
distribution using only the two scalar parameters 5 and 
82- 
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FIG. 3: Centered and reduced Probability Density Function 
of 8(t), (a) for the TM73(— ) configuration with annulus at 
seven different Reynolds numbers ranging from 1.25 x 10 s 
to 5 x 10 5 , and (b) for the TM73(+) configuration without 
annulus for 16 experiments performed with different values of 
rotation number 6 = (fi — + /a) ranging from to 

1. 8* — (S — 8)/ 62 is the centered reduced value of 8. The 
dash-dotted lines are Gaussian functions of zero mean and 
unit variance. 



A. Basic properties of 8 

1. Statistical properties 

As discussed in Sec. HJ S(t) widely fluctuates in time. 
However, its statistical properties in the exact counter- 
rotating case are rather simple since its probability den- 
sity function (PDF) is nearly Gaussian as illustrated in 
Fig. GH a )- Actually, this property holds at different 
Reynolds numbers (cf. Fig. [3Ja)) and for all forcing 



2. Dependence on the Reynolds number 

For a given forcing, we observe a very weak dependence 
of 5 with Re in the studied range as one can see in Fig. 
21 The largest dispersions, of the order of 8%, are ob- 
served in the negative rotation sense for both TM60 and 
TM73 impellers without annulus. For any other forcing 
condition, the dispersion is less than 5%. Similar conclu- 
sions hold for 82 represented in Fig. 0] using error bars. 
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FIG. 4: 5 as a function of Reynolds number Re. Negative 
Re correspond to rotation sense ( — ). (a) and (b) correspond 
respectively to TM73 and TM60 impellers, with (o), and with- 
out (□) annulus. Dash-dotted lines corresponds to the mean 
value of 8 for each configuration and 82 is represented using 
error bars. 



This behavior is not surprising since in our water ex- 
periments the Reynolds number is so large that we are 
in the fully developed turbulent regime where statistical 
properties have already been shown not to depend on 
the Reynolds number [13|, HH EH • On the contrary, we 
expect a different behavior when lowering the Reynolds 
number (l3| . since, in the limit of very low Re, the in- 
stantaneous flow is laminar and identical to the mean 
flow: 8(t) =6=1. The study of the transition be- 
tween laminar and fully turbulent regime requires the 
use of glycerol-water mixing and is beyond the scope of 
the present paper. Note however that this transition has 
already been studied in the von Karman experiment — 
TM60(-) without annulus— by Ravelet et al. Q us- 
ing a one point measurement Laser Doppler Vclocimetry 
system. It has been shown that the local value of the 
variance of the azimuthal velocity, v 2 — v 2 e , behaves as 
(Re — Rec) 1 ! 2 and saturates above Ret = 3300, where 
Re c = 330 corresponds to the transition from steady 
to oscillatory laminar flow and Ret to the onset of the 
fully turbulent incrtial regime. A similar behaviour is en- 
countered in Direct Numerical Simulations of the Taylor- 
Green flow [IH where an increase of 5 from 1 at low Re 
to a saturation value of 3 for Re above 10 3 has been 
observed. 



B. Influence of forcing and annulus 

It is difficult to estimate turbulence intensity by simply 
looking at an instantaneous velocity field (cf. Figs, [ljc) 
and (f)). To achieve this, the local turbulence intensity i 



maps (cf. Figs. [TJa) and (d)) constitutes a good qualita- 
tive tool to estimate the overall turbulence intensity and 
structure of the considered flow. However, a quantity like 
8 is needed to provide a global quantitative tool. 



1. Effect of impeller and vessel geometry 

As reported in Table [H the values of 8 and 82 allow 
to quantify the differences in the turbulence fluctuations 
level as a function of the forcing configurations (TM60 
and TM73, with or without annulus). The major result 
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0.24 
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1.52 


1.47 


1.50 


1.48 


0.11 


0.15 


0.11 


0.12 



TABLE I: Values of 8 and 82 for various configurations. 

of Table U is the effect of the annulus which systemati- 
cally reduces 8 and 82, meaning that the instantaneous 
flow is much closer to the mean flow. A great part of this 
drop probably reflects the reduction of the slow fluctua- 
tions of the shear layer already known from time spectral 
analysis to be responsible for a major part of the energy 
fluctuations [Til ]. From a spatial point of view, this re- 
duction is the consequence of the expected locking of the 
shear layer in the annulus plane which appears in direct 
visualizations of the flow. Additionally, we note that the 
annulus tends to collapse all values of 8 to 1.5 ± 0.03. In 
contrast, a dispersion of about thirty percents remains 
on the variance 82- 



2. Qualitative interpretation 

To understand the measures presented in the last para- 
graph, we report visual observations of the flow seeded 
with bubbles. Such visualizations allow to unveil spatial 
structures of turbulence and especially the largest pat- 
terns which exhibit the slowest time dynamics. Let us 
describe our observations in the TM73 impellers config- 
uration. 

Without annulus, the main structure of the free shear 
layer at very high Reynolds number consists of three big 
fluctuating vortices, i.e. a m = 3 azimuthal wave number 
mode (cf. Fig. The size of these vortices is almost 
the full free height of the vessel for sense (— ) and just a 
little bit smaller for sense (+). These vortices fluctuate 
in azimuthal position as well as in amplitude or apparent 
size. Nucleation and merging frequently occur. 

In the presence of an annulus, each vortex splits into a 
pair of smaller co-rotating vortices attached to the lead- 
ing edges of the annulus (sec the schematic view in Fig. 
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FIG. 5: Sketch of the shear-layer patterns, on the developed 
cylinder (9,z), from visual observations. Without annulus 
(top sketch), the shear produces a triplet of large corotating 
vortices, corresponding to a m = 3 mode. In the presence of 
the equatorial annulus (bottom sketch), the three vortex cores 
split into corotating pairs: there is a strong energy transfer 
from the fundamental m = 3 mode to its first m — 6 harmon- 
ics. 

[5]). The m = 6 harmonics appears and pumps an im- 
portant amount of energy to the fundamental m = 3 
mode. The strong reduction in 8 and 82 is probably the 
trace of this scale change in the shear layer structure. 
Furthermore, we observe the temporal dynamics of the 
shear layer vortices which is different for both senses of 
rotation. For sense (+), the vortices keep high mobility 
and fluctuation levels. On the contrary, for sense (— ), the 
vortices are almost steady: merging, nucleation and large 
azimuthal excursions are inhibited. The only remaining 
dynamics is an azimuthal vibration of their center at a 
frequency of a few Hertz. This qualitative observation 
is corroborated by the behavior of <5 2 : even if 8 ~ 1.5 in 
both cases, ($2 is lower in the (— ) senses, where the vortex 
pairs are more stable. 

3. Scale by scale characterization 

By construction, 8 accounts for the total fluctuating 
kinetic energy of the turbulent flow. The use of PSD 
as an intermediate step to compute 8 (cf. Eq. ([2])) can 
provide useful information about the distribution of the 
energy over the various length scales of the flow as illus- 
trated in Fig. [HI Indeed, in Fig. [Ufa), we compare the 
time-averaged PSD of the instantaneous velocity fields, 

E(k), to the PSD of the corresponding time averaged 
field, E (k), for the two TM73(+) setups, with and with- 
out annulus. 

In the spirit of the definitions of turbulent intensity 
i and of the quantity 8, and because we do access the 
full spatial spectra, we choose to normalize the PSD by 
the total kinetic energy of the mean flow. Therefore, 
the integrals of the displayed curves for the two time- 
averaged flows are equal to 1. The wave number k is 
given by: k 2 = k 2 r + k 2 . 

In the time-averaged PSD of the instantaneous flow, 
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FIG. 6: (a) Two-dimensional time-averaged Power Spec- 
tral Density of instantaneous flows (solid lines) and Power 
Spectral Density of the corresponding time-averaged flows 
(dash-dotted lines) for two experiments performed with setup 
TM73(+), without the annulus (□) and with the annulus (o). 
The PSD are normalized by the total kinetic energy of the 

mean flow, (b) 8(k) as function of k for the same setup with- 
out (□) and with annulus (o). The specific data point 6(1/ D) 

is represented with a black symbol (■). The difference A8(k) 
between the two curves is also displayed (plusses). At high 

k, it is of the order of 5(1/ D) without annulus (see text for 
details). Wavenumbers are normalized by the PIV window 
size D = 1.93.R. 



the major part of the fluctuating kinetic energy is con- 
centrated close to the forcing scale, where the energy is 
injected, like in any homogeneous and isotropic turbu- 
lent flow. On the contrary, at higher k, the slopes of the 
different PSD are close to —2 and larger than the clas- 
sical —5/3 exponent observed for homogeneous isotropic 
turbulence. 

At first sight, the time-averaging removes kinetic en- 
ergy at all spatial scales. A closer look indicates different 
relative reductions depending on the presence or the ab- 
sence of the annulus. At high k, the ratio of the PSD of 
instantaneous to time-averaged flows is of the same order 
of magnitude for both configuration. However, at low fc, 
this ratio is 2 (resp. 1.2) without (resp. with) annulus: 
there are five time less fluctuating low-spatial-frequcncics 
to be averaged when the annulus is inserted meaning that 
the instantaneous flow is closer to the time-averaged flow. 
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To quantify more precisely the weight of each scale in 
S(t), we introduce S(k,t) defined as: 



8(k,t) 



£dk'E(k',t) 
J~dk'E (k')' 



so that S(k,t)—>S(t) when k — > oo. In Fig. EJb), S(k) 
is plotted for TM73(+) flows with and without annulus. 
These curves represent the integral from to k of the 
above normalized PSD and they converge towards 8 at 
large k. 

First of all, we see that the contribution of the 5 first 
wave numbers, i.e. of spatial modes larger than D/5 ~ 
R/10, is about eighty percents of the total value of 8: 
again, large scales dominate the fluctuations. 

Some quantitative analysis can be carried on. We note 
that the fluctuation level in the first k = 1/D mode of 

the 8(k) curve without annulus and the difference A<5 be- 
tween the two flows — with and without annulus — are of 
the same order of magnitude. This is visible in Fig.[6jb), 

where the 6(1/ D) data point has been emphasized, and 
has been verified for the three other couples of flows. As 
described in the previous section, the effect of the annulus 
can be seen as a spatial filtering of the largest patterns of 
the flow. The difference between the two curves can be 
described by a simple empirical protocol. The amount 

of 6 in excess, i.e. 5(1/ D) ~ A<5, is subtracted from the 

curve S(k) without annulus in a geometric progression 

over k, e.g. S(l/D)/2 k for each wavenumber k starting 
from D~ x or 2D -1 depending on the configuration. The 

result is a curve really close to S(k) with annulus. The 
effect of the annulus may thus be seen as a high-pass 
filtering, which confirms the sketch presented in Fig. [5] 
Finally, the quantity 8(k,t) appears as an efficient tool 
to track the spectral changes in turbulent scales. 



lation cells when \9\ < 9 C . A good quantity that al- 
lows to localize this position is the zero iso-surface of the 
stream function ip(r, z) of the mean flow defined through 
0,mJ) = V x (r ipeg) in cylindrical coordinates. 
Thus, in Figs. [7ja) and (b), we have plotted two mea- 
sures of the shear layer vertical position at r = 0, i.e. 
the stagnation point, and at r = 0.7 on the side. The 
discrepancies between these two sets of data renders the 
fact the shear layer is in general a curved surface expect 
for the 9 = configuration. Actually, its vertical position 
at large r is always closer to the z = equatorial plane of 
the vessel than at r = 0. As \6\ increases from zero, the 
shear layer is attracted by the slowest impeller. With- 
out annulus (cf. Fig. EJa)), it results in a global drift of 
the shear layer [TtI ] accompanied by the appearance of a 
moderate curvature. With annulus (cf. Fig. [7[b)), the 
outer edge of the shear layer remains pinned on the an- 
nulus whereas its center — the stagnation point — drifts, 
increasing the layer curvature up to a high level just be- 
low 9 C : the annulus strongly stabilizes the shear layer 
near the equatorial plane and maintains it to higher ro- 
tation numbers. 

At 9 = 9 C , we observe a discontinuity of the shear 
layer vertical position which corresponds to the turbu- 
lent bifurcation. The discontinuous bifurcation does not 
present, in the two studied cases, any hysteresis in 9. 
However, close to 9 C , we can observe transitions between 
two metastable one cell and two cells states which follow 
a slow dynamics [28| . 

In this section, we first analyze how the parameters 
d and 82 behave with the rotation number. Thereafter, 
using these tools computed over the two symmetric half 
of the flow, we analyze how they can also provide a proper 
characterization of the symmetry breaking. 



1. Variation of 6 with the rotation number 



C. Properties of 8 when /1 7^ fo 

As already mentioned in Sec. IIII A 3[ depending on the 
value of the rotation number 9 = (fi — f%)/(fx + f-z), the 
mean von Karman flow exhibits two different topologies. 
Indeed, for \9\ < 9 C , the mean flow is composed of two 
toric recirculation cells separated by an azimuthal shear 
layer, as for \9\ > 9 C , it is composed of a single recircula- 
tion cell Q . We have performed a set of experiments for 
the specific forcing condition TM73(+) with and with- 
out annulus, varying 9 from —1 to 1. In this geometry, 
9 C = 0.17±0.01 with annulus and 9 C ~ 0.095±0.005 with- 
out annulus. The aim of these experiments is to study 
the influence of the global rotation and of the transition 
between the two flow topologies on turbulence intensity. 

A way to quantify precisely the turbulent bifurcation 
of the mean flow occurring at 9 C is to study the position 
of the azimuthal shear layer that separates the recircu- 



The parameter 8 can be used to study the changes in 
turbulence intensity related to the turbulent bifurcation 
undertaken by von Karman flows at critical 9 C . Figs. [7| 
(c-f ) show the variations of 8 and 82 , calculated over the 
whole flow function 9. 

First of all, for \9\ < 9 C , 8 and 82 are relatively high 
(8 ~ 2 and 82 — 0.2). This reflects the presence of the 
highly fluctuating shear layer in the flow bulk. On the 
contrary, for \9\ > 9 C , i.e. for a large rotation-to-shear 
rate, 8 and 82 are quite small (8 ~ 1.2 and 82 — 0.05) 
and are almost independent of 9 [2{|. This means that 
the level of global rotation has no significant influence on 
the turbulence intensity level unless the mean flow topol- 
ogy changes. Actually, 8 and 82 carefully trace back the 
flow topology changes induced through the bifurcation 
at ±9 C . They might even be used as order parameters 
of such turbulent bifurcation and they provide a reliable 
measurement of threshold 9 C . 
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FIG. 7: (a) and (b), vertical position of the shear layer at r — on the rotation axis (o) and at position r — 0.7 (•), (c) and (d) 
S, and (e) and (f) 82, as a function of the rotation number 9. Left and right columns correspond respectively to setup TM73(+) 
without and with annulus. In (a) and (b), the (o) markers correspond to the stagnation point vertical position. In (c-f), 8 and 
82 are computed over the whole flow. 



2. Effect of the annulus 



Quantifying the symmetry breaking 



We have seen how the annulus stabilizes the shear layer 
and moves the bifurcation threshold. It has also a strong 
effect on the turbulence level. In Figs. [7] (c-f), we ob- 
serve that for \9\ < 6 C , when the flow is composed of 
two counter-rotating toroidal cells, the turbulence level 
is much larger without than with the annulus: the re- 
sults for zero-rotation-number (cf. Sec. IIVB[) extend to 
-8 C <6< 6 C . 

We have already mentioned that this strong reduction 
in 8 and 82 with the annulus is due to the stabilization, es- 
pecially at large scales, of the shear layer that is trapped 
by the annulus. Now, picturing that the flow has a given 
energy gap to overcome in order to bifurcate from two 
cells to one cell, the lower level of fluctuations in the 
case with annulus explains why it is necessary to explore 
larger 9 before turbulent bifurcation occurs, and then, 
why 9 C is larger with the annulus. Actually, the annulus 
is postponing the bifurcation in 9. 

For \9\ > 9 C , the low turbulence intensity, almost un- 
changed at first order, is slightly larger in the presence 
of the annulus: the one-cell flow is almost insensitive to 
the annulus presence and the slight increase is probably 
due to the vertical step flow over the annulus. 



For further exploration of the turbulent bifurcation 
and related symmetry breaking, we now calculate 5 and 
82 over half vessels, i.e. on each side of the annulus equa- 
torial plane (cf. Fig. [5]). 



Starting from the exact counter-rotating configuration, 
as \9\ is increased from zero, 8 and 82 are becoming larger 
for the half part of the flow corresponding to the slowest 
impeller and are decreasing for the other half part. This 
is directly connected with the position of the shear layer 
which gets closer and closer to the slowest impeller when 
9 increases (cf. Figs. [Tfa) and (b)). This displacement 
goes on until the level of fluctuations is strong enough to 
activate the transition of the flow from two to one cell at 



Indeed, we clearly see in Figs. [5] the peaks in 8 and #2 
that correspond exactly to the bifurcation point 9 = 9 C . 
With the annulus, where the bifurcation is much sharper, 
one could detect a critical divergence for (8 — 1) (cf. Fig. 
!>)). 
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FIG. 8: (a) and (b) S, and (c) and (d) 82, as a function of the rotation number 6. Left and right columns correspond respectively 
to setup TM73(+) without and with annulus. (a) correspond to 8 computed over the upper half part of the flow as (V) to the 
lower half part. 



V. CONCLUSIVE DISCUSSION 

We have introduced a new global quantity 8 that 
characterizes turbulent fluctuations in inhomogcneous 
anisotropic flows. This time-dependent quantity is based 
on spatial averages of global velocity fields rather than 
classical temporal averages of local velocities. We have 
shown that, generally, properties of 8(t) seems to be fully 
provided thanks to only two parameters, its time aver- 
age 8 and variance 82- These parameters generalize to 
inhomogcneous anisotropic flows the classical notion of 
turbulence intensity, based on local, single point mea- 
surements. 

Properties of 8 and 82 have been experimentally stud- 
ied in the typical case of the von Karman flow for differ- 
ent forcing and geometries. We have shown that, in the 
fully turbulent regime, they arc Reynolds independent, 
like any classical inertial range quantity. However, 8 and 
82 depend on forcing and geometry and faithfully reflect 
major changes in the flow topology. They can therefore 
be used as a new tool for comparison of different turbu- 
lent flow. In the present paper, we provided an example 
in which 8 is used to characterize the turbulent bifurca- 
tion in the von Karman flow induced through differential 
rotation of the two impellers. Finally, 8 and 82 are used 
to characterize the scale by scale energy budget as a func- 
tion of the forcing mode as well as the transition between 
the two flow topologies. 

Another interesting application would be the study of 



the generation of magnetic field by a turbulent flow, the 
so-called dynamo instability. This problem has attracted 
recently a lot of experimental attention [l8|, [l9|, [2(| ■ In 
the case where the corresponding turbulent flow has a 
non-zero mean value, the dynamo instability may be seen 
as a classical instability problem, governed by the mean 
flow and the fluctuations [l|| [2l|, [24] . A natural question 
in this case is therefore how to quantify the relative level 
of fluctuations, and the deviations from the mean flow 
they induce, so as to implement efficient control strate- 
gies to decrease or increase their influence. Since our 
global parameter quantify the difference between the in- 
stantaneous and the mean flow, they are natural can- 
didate to discriminate between different dynamos. This 
has been recently illustrated in numerical simulations of 
Taylor-Green flow (23|. No equivalent quantitative analy- 
sis has been performed for experimental dynamos. How- 
ever, among these successful experimental dynamos, we 
can suspect that the Karlsruhe and Riga [l8[ flows 
are characterized by values of 8 lower than in von Karman 
flows. In the future, we plan to use these parameters for 
the analysis of recent results of the von Karman Sodium 
(VKS) experiment. 



APPENDIX A: CONVERGENCE TOWARDS THE 
MEAN FLOW 

As we have seen in Sec. [Hi the parameter 8 can be seen 
as the average square distance between the instantaneous 
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FIG. 9: 8—1 computed with velocity fields averaged over 
A^ instantaneous fields for various setups : TM73( — ), with- 
out the annulus (■) and with the annulus (•) and TM73(+), 
without annulus (□) and with annulus (o). A guide to the eye 
for the N~ J dependency is plotted with a dotted line. 

and the time averaged velocity fields. Therefore, we can 
use a slight modification of 5 to study the convergence 
towards the mean flow through statistical averaging. For 
this, we define V(N, t) the velocity field averaged around 



time t over N instantaneous fields. N < N max where 
N m ax = 5000 is the total number of instantaneous fields. 
From a practical point of view, V{N max ) = V(t). Then, 
we define: 



5(N) = < y W». (Al) 
(V ) 

With this definition, 5(1) = 5 and S(N max ) = 1. More- 
over, S(N) — 1 measures the average square distance be- 
tween the partially averaged field V(N) and the mean 
flow V, so that its variations with N can be used to 
study the convergence towards the mean flow. We can 
see in Fig.|H]that this square distance, S(N) — 1, decreases 
as A^ _1 , at least at large N, what is typical of an uncor- 
rected fluctuating quantity. However, for the negative- 
rotation-sense-with-annulus case, the A^ -1 dependency 
is observed only at N larger than 100. Actually, in that 
particular setup, we have observed, by means of bubble 
seeding (cf. end of Sec. IIVB 2p . that even if the largest 
structures of the shear layer where removed by the annu- 
lus, a pair of coupled vortices appeared. These vortices 
are smaller than those we observe without the annulus, 
but their long time coherent precessing must induce long 
time correlations that slow down the decrease of 6(N) . 
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